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Abstract 

Dispersion relations allow for a coherent description of the nucleon electro- 
magnetic form factors measured over a large range of momentum transfer, 
Q 2 ~ . . . 35 GeV 2 . Including constraints from unitarity and perturbative 
QCD, we present a novel parametrisation of the absorptive parts of the var- 
ious isoscalar and isovector nucleon form factors. Using the current world 
data, we obtain results for the electromagnetic form factors, nucleon radii 
and meson couplings. We stress the importance of measurements at large 
momentum transfer to test the predictions of perturbative QCD. 



# 1 Work supported in part by the Deutsche Forschungsgemeinschaft (SFB 201). 



I. INTRODUCTION AND SUMMARY 



The electromagnetic structure of the nucleon as revealed in elastic electron-nucleon scat- 
tering is parametrized in terms of the four form factors Fi£(Q 2 ) (with Q 2 the squared mo- 
mentum transfer). The understanding of these form factors is of utmost importance in any 
theory or model of the strong interactions. Abundant data on these form factors over a 
large range of momentum transfer already exist, and this data base will considerably im- 
prove in the few GeV region when experiments at CEBAF will be completed and analysed. 
In addition, experiments involving polarized beams and/or targets are also performed at 
lower energies to give better data in particular for the electric form factor of the neutron, 
but also the magnetic proton and neutron ones. Such kind of experiments are under way 
at MAMI, ELSA, MIT-Bates and other places. Clearly, theory has to provide a tool to 
interpret these data in a model-independent fashion. Many years ago, dispersion theory 
was developed to extract these form factors from elastic ep and ed scattering data (and 
others) (for some early references, see e.g. @ ||.). Such a dispersion theoretical analysis 
is largely model-independent with the exception of the absorptive parts of the form factors, 
which are often parameterized in terms of a few meson poles. A large class of the early 
models of the spectral functions related to the imaginary parts of Ff ^{Q 2 ) wa s based on 
the successful vector meson dominance (VMD) hypothesis (for reviews, see [Q ||) which 
states that a photon couples to hadrons only via intermediate vector mesons. However, as 
already pointed out in 1959 by Frazer and Fulco || 0, such an approach is at odds with 
general constraints from unitarity. The two-pion continuum has a pronounced effect on the 
isovector spectral functions on the left wing of the p-resonance. This becomes particularly 
visible in the determinations of the corresponding nucleon radii. This effect was quantified 
by Holder and Pietarinen |§ but has been neglected in most of the recent discussions of the 



nucleons form factors, like in M [10| [11| [12| [ 13| . For an accurate determination of the 



isovector nucleon radii, the inclusion of the two-pion cut contribution is mandatory. An- 
other constraint comes from the accurate measurement of the neutron charge radius in very 
low-energy neutron-atom scattering (for a recent update, see e.g. Leeb and Teichmeister 
|14|). This leads to an additional normalization condition thus reducing the number of free 
parameters.^ Furthermore, in the framework of perturbative QCD (pQCD), the behaviour 
of the pertinent form factors for large momentum transfer can be inferred up to a model- 
dependent normalisation [JTE| . This has lead to models which try to synthezise aspects of 
low-energy hadron physics with the ones from pQCD, see e.g. |[10| [|ITJ |TjJ |18| fSOfl . No 



clear picture has yet emerged at which momentum transfer the asymptotic behaviour really 
sets in. We stress here that any serious attempt to describe the electromagnetic structure 
of the nucleon has to account for all these constraints. 



# 2 One can also consider the determination of the electric charge radius of the proton by Simon 
et al. [15] as a further low-energy constraint. It is, however, less stringent as the one related to 
the neutron charge radius and will be discussed separately in section II E . 
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The work presented here is an update and extension of the classical paper by Holder and 
collaborators ||21|| . As in that paper, we determine the p-meson contribution to the isovector 
form factors from extended unitarity as discussed by Frazer and Fulco |7j . While we use the 



old iiN partial waves extrapolated to the time-like regime 22], we update the pion form 



factor to account for the measured p — uo mixing |[23|| . In addition, since we are dealing 
with what is called an ill-posed problem, we restrict the number of additional pole terms 
in the isovector and isoscalar channel to a minimum (see the discussion below). We also 
implement the large-Q 2 behaviour dictated by pQCD, similar in spirit to the work of Furuichi 
and Watanabe [T7J |TJ|. However, as we demonstrate later on, their choice of the spectral 
functions is not compatible with pQCD. What one ends up with is a set of superconvergence 
relations which reduce the number of free parameters. These superconvergence relations can 
be implemented in a variety of ways. Therefore, we will discuss two versions of the spectral 
functions with the proper asymptotic behaviour guided by simplicity and the notion of 
separating the hadronic from the pQCD contribution. A fit to the accessible data then 
allows not only to extract the nucleon electromagnetic form factors, but also the nucleon 
radii and nucleon-meson coupling constants. In particular, we also find a large coupling of 
the $-meson to the nucleon, seemingly in contradiction with the OZI-rule as stressed in 



the work by Genz and Holder [2l] and, more recently, by Jaffe [25]. It should be already 



noted here that this large $-nucleon coupling can be understood in terms of coupling of the 



photon to the kaon cloud surrounding the nucleon [26]. Our approach furthermore sheds 
light on the possible onset of pQCD and the regime where data are urgently called for. For 
a status review on the determination of the nucleon electromagnetic form factors, we refer 
to Holder p], Milner P§ and Bosted 



The pertinent results of this investigation can be summarized as follows: 

(i) Including all the constraints discussed and using the current world data, we have found 
a new fit to the nucleon electromagnetic form factors. Besides the two-pion continuum 
and the pQCD contributions, we have three isoscalar and three isovector poles. This 
is the minimum number required by the data. We also stress that it is mandatory to 
simultaneously fit the proton and the neutron data. 

(ii) The form factors F[ s \t) and F^\t) show a marked dipole structure. The reasons is 
that the two lowest poles in the respective channels are not too far separated in mass 
and appear with residua of the same magnitude but different in sign. In the isoscalar 



case, this is due to the closeness of the u and $ mesons [21 1. The novel feature of 



our study is that the first pole above the two-pion continuum can be identified with 
the //(1450) leading to the dipole structure in (t). We furthermore find the third 
isoscalar and the second isovector pole at M s > = 1.60 GeV and M p u = 1.65 GeV. Only 
the third isovector pole cannot be identified with a physical state. 
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(iii) We have given a new determination of the electric and magnetic radii of the proton 
and the neutron. We find r p E = 0.847 fm, r p M = 0.853 fm and r\ t = 0.889 fm, all with 
an uncertainty of about 1%. These are consistent with the findings of ref. with the 
exception of the neutron magnetic radius. The difference is mainly due to the neglect 
of the super convergence relation related to the asymptotic behaviour of in [21|. 



(iv) We have shown that the accurate determination of the proton charge radius in ref. [15 
is only consistent with the dispersive analysis within one standard deviation. Using 
the central value of ref. |15| for r E , one cannot simultaneously fit the proton and the 



neutron data. 

(v) We have determined the vector (gi) and tensor (g 2 ) couplings of the to and the $ 
mesons. The gi are slightly larger than in Ref. EIJ while the g 2 are of comparable 
size. We stress again that the large &NN coupling does not indicate a violation of the 
OZI-rule but rather accounts for the neglect of the KK continuum. 

(vi) In our best fit, the parameter A 2 , which is a measure of the boundary between the 
hadronic (pole) and quark (pQCD) contributions, comes out to be A 2 ~ 10 GeV 2 . 
Only in the case of the proton magnetic form factor one has data at suffciently 
large momentum transfer to possibly test the predictions of pQCD. Our conclusion 
is that even in that case more data are needed to really see the asymptotic be- 
haviour Q 4 G P M (Q 2 ) — > constant (modulo logarithms). For all form factors there 
is still a sizeable hadronic contribution in the momentum transfer range between 
Q 2 = 10 ... 20 GeV 2 . In particular, the data available for the ratio Q 2 F 2 {Q 2 ) / F^Q 2 ) 
are at too low momentum transfer to indicate any scaling. 

(vii) The fits are completely insensitive to the number of flavors entering via the anoma- 
lous dimension appearing in the asymptotic behaviour of the form factors and to the 
p — uj mixing which comes in via the calculation of the isovector two-pion unitarity 
correction. 

(viii) More accurate data at low, intermediate and (very) large momentum transfer are called 
for to further tighten the constraints on the radii and coupling constants and to really 
test the pQCD predictions. 

The manuscript is organized as follows. In section |I| we discuss the basic concepts 
underlying the dispersion-theoretical analysis of the nucleon electromagnetic form factors. 
In particular, we summarize the various constraints which have to be implemented to obtain 
a consistent picture. These are the inclusion of the two-pion continuum in the isovector 
spectral functions (sect. |II C| ), the determination of the neutron charge radius from very low 
energy neutron-atom scattering (sect. [II D|) and the behaviour of the various form factors as 



given by pQCD (sect. [11 F| ). We also discuss another less stringent constraint related to the 
proton charge radius (sect. [II E| ). Section [HI] contains the form factor parametrizations we 



will use. We present two forms, one in which the pQCD constraints are implemented in the 
most simple fashion and another one, which allows for a separation between the hadronic 
(pole) and the pQCD terms already on the level of the spectral functions. The results are 
presented and discussed in section |V[ The appendices A and B contain some technicalities, 
and in appendix C, we give the explicit parametrization of our best fit. 



3 



II. BASIC CONCEPTS 



In this section, we assemble all tools necessary for the construction of the fit functions for 
the spectral distributions discussed in section |T|. We also supply the basic definitions and 
notations underlying our analysis. In particular, we discuss the various low and high energy 
constraints which have to be imposed and allow to reduce the number of free parameters. 
This material is mostly not new but necessary to keep the manuscript self-contained. The 
reader familiar with it is invited to skip this section. 



A. Definition of the nucleon electromagnetic form factors 

The matrix-element of the electromagnetic (em) current sandwiched between nucleon 
states can be expressed in terms of two form factors, 

< N(p') | J, | N(p) >= e 1 , F 1 (t) + ^— a, v Q v F a (*)} u(p) , (1) 

with rriN the nucleon mass, t = (p' — p) 2 = —Q 2 the invariant momentum transfer squared 
and t < in the space-like region. Fi(t) and F 2 (t) are called the Dirac and the Pauli form 
factor, respectively. They are normalised at t = as 

Ff(0) = l, Ff(0) = 0, F 2 p (0) = Kp , F 2 "(0) = K n , (2) 

with K p (n n ) the anomalous magnetic moment of the proton (neutron) in units of the nuclear 
magneton, [i N = e/2rriN- For later use, we also need the isospin decomposition, i.e. the 
response of the nucleon to the isoscalar (denoted by the superscript '( s )') an d isovector 
(denoted by the superscript \v) J ) components of the electromagnetic current, 

^ (S) = \{Ff + F?) , Ff> = \(Ff - F?) , (2 = 1,2), (3) 
with the normalisation 



F/(0)=z/, i = l,2; /=(«),(«), (4) 

where 

U [ S) = = 1 , U { 2 S) = ]-{n p + K n ) , = l(Kp - «n) . (5) 

One can also introduce the so-called Sachs form factors G p ^ n M (Q 2 ) related to the F\{Q 2 ) and 
F 2 (Q 2 ) via 

G e (Q 2 )=F 1 (Q 2 )-tF 2 (Q 2 ), G m (Q 2 )=F 1 (Q 2 )+F 2 (Q 2 ), t = -^-. (6) 
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In the Breit frame, Ge and Gm are nothing but the Fourier-transforms of the charge and 
the magnetization distribution, respectively. 

The slope of the form factors at t — is conventionally expressed in terms of a nucleon 
radius < r 2 >^ 2 , 



Fit) = /•'(()) ( 1 + g < r 2 > t + 



(7) 



which is rooted in the non-relativistic description of the scattering process, in which a point- 
like charged particle interacts with a given charge distribution p(r). The mean square radius 
of this charge distribution is given by 



< r >- 



dr 4nr 2 p(r) 



6 dF{Q 2 



F(0) dQ 1 



Q 2 =0 



(8) 



Eq.(||) can be used for all form factors except G% and F™ which vanish at t = 0. In these 
cases, one simply drops the normalization factor 1/F(0) and defines e.g. the neutron charge 
radius via 



< r 2 >? 



-6 



dG E (Q 2 ) 



dQ 2 



(9) 



To conclude this section, we remark that the slopes of G n E and F™ are related via 

dG E (Q 2 ) dF?(Q 2 ) F?(0) 



dQ 2 



Q 2 =0 



dQ 1 



4:171% 



(10) 



The second term in eq. (|T0|) is called the Foldy term. It gives the dominant contribution to 
the slope of G E . 



B. Dispersion relations for the nucleon em form factors 

Let F(t) be a generic symbol for any one of the four nucleon em form factors Ffy if). 
We assume the validity of an unsubtracted dispersion relation of the formp 3 ! 



F(t) 



ImF(t') 

7T J to f — t 



dt' 



(11) 



The spectral function Im F(t) is different from zero along the cut from to to oo with to = 
4 (9) M 2 for the isoscalar (isovector) case and denotes the pion mass. The proof of the 
validity of such dispersion relations in QCD has not yet been given |30|. Eq. (|i~ID means 
that the em structure of the nucleon is entirely determined from its absorptive behaviour. 



* 3 Since the normalization F(0) is known, one could also work with once-subtracted dispersion 
relation. 
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Data for F(t) are given for t < 0. If these would be infinitively precise, the continuation to 
other values of t > 4M 2 . would be unique. However, the available data have certain error 
bars which make the continuation unstable as detailed in the review |31]]. This is what one 



calls an ill-posed problem. Therefore, one needs some additional assumptions which should 



have a physics motivation. We follow here the prescription due to Sabba-Stefanescu |32 
and require that the spectral functions should not have more than the minimum number of 
oscillations and parameters required by the data. To make this more specific, let us consider 
as an example the rather successful dipole fit, 

G p E {t) ~ fx P G p M (t) ~ imG n M {t) ~ G D (t) , G D (t) = (1 - t/0.71 GeV 2 )- 2 , (12) 

where fj, p and \i n are the magnetic moments of the proton and the neutron, respectively. 
Such a form is only possible if the corresponding spectral functions contain at least two pole 
terms with opposite signs. 

Of particular interest is the VMD approach in which the spectral functions are approx- 
imated by a few vector meson poles, namely the p, . . . in the isovector and the u, $, . . . in 
the isoscalar channel, in order. In that case, the form factors take the form 



v lvl (i) 



F * I) ( t ) = Y,J^ L — t i = l,2; I = v,s, (13) 



with 



/If 2 

V _ IVI V n VNN 



a;=^-gr N , V = />M...,i = l,2. (14) 



fv 

In many fits, the mass parameters My are taken from the known meson masses and the 
strength parameters a( are fitted. This in turn leads to a determination of the various 
vector-meson-nucleon coupling constants, g( NN ■ Note, however, that such a procedure 
becomes increasingly arbitrary for higher mass excitations. One defines the ratio of the 
tensor to vector coupling, Ky, via 

gVNN 

K V = n VNN ■ ( 15 ) 

yi 

One expects e.g. k p to be large (~ 6) and to be small (~ 0). In strict VMD, one has 
k p = k p — K n = 3.71 and = k v + n n = —0.12. The coupling constants fy are known from 
the widths of the leptonic decays V — > e + e~, i.e. 



fy « 2 M v 



(16) 



4tt 3 T(V -> e+e~) 
Clearly, such pole terms contribute to the spectral functions in terms of 5-functions, 

lmF^(t) =iraj 5{t-M v ) . (17) 

Of course, there are other contributions related to intermediate states like mr (n > 2), NN, 



KK and so on. As we will discuss in section [II Q , of these the 2tt intermediate states play 
the most important role. 
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C. Constraints from unitarity 



Here, we briefly summarize what is known about the contribution of the two-pion contin- 
uum to the isovector spectral functions and how this should be implemented. The unitarity 
relation of Frazer and Fulco determines the isovector spectral functions from t = AM 2 to 
t ~ BOM 2 ~ 1 GeV 2 in terms of the pion form factor F n (t) and the P-wave inrNN partial 
wave amplitudes, cf. Fig.l. We use here the form 

lmFt ) ' 2n (t) = ^=\F T (t)\ 2 J l (t), i = l,2 (18) 

with q t = Jt/4: — M% and the functions Ji(t) are related to the P-wave nN partial waves in 
the t-channel, f±(t), via 

fl(t)=F n (t)J ± (t) , (19) 



in the conventional isospin decomposition. The Jj(i) are tabulated in [Q. For the pion form 
factor, we use the recent parametrization of Barkow et al. |£J which takes into account p — u 
mixing, 

m=K{t) l \ a f" {t) (20) 

1 + a w 

with a w the mixing parameter, a u = 0.0038 [^3|]. The functions F^(t) (V = p,u) are of the 
standard Gounaris-Sakurai form |Q . A typical result for the corresponding isovector spec- 
tral functions (weighted with 1/t 2 ) is shown in Fig.2. One notices the strong enhancement 
close to two-pion threshold. The reason for this behaviour has been known for long. The 
partial waves f±(t) have a branch point singularity on the second sheet (from the projection 
of the nucleon pole terms) located at 

t c = AMI ~ M 4 Jm 2 N = 3.98 M 2 n , (21) 

very close to the physical threshold at to = 4M 2 . The isovector form factors inherit this 
singularity and the closeness to the physical threshold leads to the pronounced enhancement. 
Note that in the VMD approach this spectral function is given by a 5-function peak at 
M 2 ~ 30M 2 and thus the isovector radii are strongly underestimated if one neglects the 
unitarity correction || as can be seen from the formula 

<r 2 >^=H2^ ^Imi^V). (22) 

Consequently, in the isovector channels one should not use a simple p pole but rather work 
with the spectral functions as given by eq. ([Op . This is the procedure adopted here. To 
speed up the numerical calculations, one can fit the corresponding two-pion contributions 
to the isovector form factors via (as it was done in plf) 
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F^{t) = Fnt) = a ^^:^~ 2 '\ i = l,2 (23) 

with a{ = 1.0317, a p 2 = 5.7824, b{ = 0.0875, = 0.3907, c{ = 0.3176, c£ = 0.1422, 
d{ = 0.5496 and d p 2 = 0.5362. We have also performed fits with the exact representation 



eq.([18D which lead to the same results as the use of the form eg . fl23|) . To end this section, we 
remark that the form eq.([]l|) is exact below four-pion threshold, t = 16M 2 and we stop at 
^up = 50M 2 since the pion form factor shows some structure on the right wing of the p-peak 
which can not simply be fitted by a superposition of Gounaris-Sakurai functions. In the 
isoscalar channel, it is believed that the pertinent spectral functions rise smoothly from the 
three-pion threshold to the u peak, i.e. that there is no pronounced effect from the three- 
pion cut on the left wing of the cu-resonance (which also has a much smaller width than the 
p). Chiral perturbation theory |35j [36] or an investigation of the spectral functions related 
to the process NN — > 3n could be used to settle this issue. The one loop calculation of the 
isovector nucleon form factors indeed shows the unitarity correction on the left wing of the p 
1 35 1 - For the isoscalar form factors, a two loop calculation of the pertinent imaginary parts 
would reveal whether there is some enhancement around t = 9M 2 or justify the common 
assumption that one has a smooth isoscalar spectral functions driven by the u at low t. 



D. Constraints from low energy neutron atom scattering 



There exists a large number of experiments trying to determine the neutron-electron 
scattering length b ne from low-energy neutron-atom scattering. The interest in this quantity 
stems from its direct relation to the mean square radius of the neutron, 

am. 



< r 2 >? 



(24) 



or 



dG 7 ^ 



(Q 2 



dQ 2 



-14.39 [fm] b r . 



(25) 



with a = e 2 /4ir = 1/137.036 the fine structure constant. For an early review, see e.g. 



Foldy |37j . There has been some controversy about the actual value of b ne over the years as 
discussed in [jTI| . We use here the most recent value, b ne = (—1.308 ± 0.05) • 10~ 3 fm 
which leads to (adding systematic and statistical errors in quadrature) 



38 



-6 



dF?(Q 2 



dQ 2 



Q 2 =0 



(0.0136 ± 0.0043) fm 2 



(26) 



using eq. ([[I]). Notice that due to the large cancellation between the Foldy term and the 
slope of G% the slope of F" is very sensitive to the value of b ne . If one uses e.g. the 
older value of Koester et al. f3j, b ne = (-1.32 ± 0.04) • 10~ 3 fm, one finds -QdF^/dQ 2 = 
(0.0126 ± 0.0035) fm 2 . The slope of F™(Q 2 ) as given in eq.(]26f) will be imposed on all our 
fits. 
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E. Constraints from low energy electron proton scattering 



Simon et al. I5| have presented a precise measurement and analysis of their and other 



existing data for elastic electron-proton scattering in the range Q 2 = 0.005 . . . 0.055 GeV 2 . 
They have performed a fit of the type 

G p E (Q 2 )=a (l+Q 2 - + Q 4 -) , (27) 



do do, 

with adjustable parameters ao, ai, Q2- The proton charge radius is thus given by 



dG p E (Q 2 



G^(0) dQ< 



.2 



a 6ao 



For the data in the low-energy region, the contribution of the Q A term to the proton electric 
form factor is marginal (< 0.3%). This leads to an rather accurate value for < r E > p , 



< r 2 E > p = (0.862 ± 0.012) fm 2 . (29) 

With that constraint, the authors of [15] performed a four pole fit (with two masses fixed 
at M p = 0.765 GeV and Mp> = 1.31 GeV) to the available data for the proton electric and 
magnetic form factors up to Q 2 ~ 5 GeV 2 . This allowed to reconstruct the spectral function 
lmG v E (t). It is in agreement with the one of ref. [§j with the exception of the region around 
t = 7M 2 . One could now argue that the dispersive analysis to be performed should be 
constrained by the value of the proton radius, eq. (|29|) . In that case one would, of course, 
have to exclude the low-energy data (Q 2 < 0.6 GeV 2 ) from the data basis. So as not to 
exclude the possibility of new and more precise data also for G P E at low Q 2 , we will perform 
a set of fits with all data included and another set with the constraint eq. (|29|) imposed and 
the corresponding data removed. We consider this as a good measure of the accuracy in 
determining the various nucleon radii from a dispersion-theoretical analysis. 



F. Constraints from perturbative QCD 

Perturbative QCD allows to constrain the behaviour of the nucleon em form factors for 
large momentum transfer, Q 2 = — t — > oo. In its most simple fashion, the so-called quark 
counting rules |4(J give the leading power in the large-Q 2 fall-off of the form factors by 



counting the number of gluon exchanges which are necessary to distribute the large photon 
momentum equally to all constituents. In the limit — t — >• oo this leads to 

{-t) i+1 Fi(t) -> constant , i = 1,2 (30) 

which in turn translates into a set of superconvergence relations for the spectral functions, 

1 r°° 

- / dt'lmF^t') = , (31) 

71 Jto 
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1 POO 1 POO 

- / dt'lmF 2 (t') = - / dt' t' lmF 2 (t') = , 

7T Jt 7T Jt 



(32) 



for both the proton and the neutron. These arguments have been sharpened in fllql . There, 
it was shown that the magnetic Sachs form factor takes the form 



Gm(Q 2 



Q 4 



lni Ql 



2 



(33) 



with a s (Q 2 ) the running strong coupling constant, 



a s (Q 2 



Air 



P HQVQl) 



(3=11 



-N f 



(34) 



Here, Nf is the number of flavors, (3 the QCD /3-function to one loop and Qo ~ Aqcd- The 
constant C in eq.(|33D is model-dependent and subject of much controversy. Its explicit form 
is not needed here. Taking furthermore into account that any helicity flip leads to an extra 
1/Q factor, one finds 



-(i+l) 



In 



7 = 2 + 



3(3 



1,2 . 



(35) 



The anomalous dimension 7 depends weakly on the number of flavors, 7 = 2.148, 2.160, 
2.173 for Nf = 3, 4, 5, in order. In what follows, we will construct spectral functions which 
lead exactly to the large-t behaviour as given in eq . (|35D . 



III. PARAMETRIZATION OF THE FORM FACTORS 

We have now assembled all tools to construct parametrizations of the nucleon form 
factors subject to the constraints discussed in the previous section. We will present two such 
parametrizations. The first one, which we label 'multiplicative', is guided by simplicity. For 
the second one, called 'additive', the spectral function is split in a way which allows for a 
clearer separation between the hadronic (pole) and the pQCD contributions. 

A. Multiplicative parametrization 

Consider the functions 



Ft\t) = Fr(t)L(t) 



E 



In 



A 2 -t 



(36) 



Ft\t) = Fr(t)L(t) 



0), 



a { i ] L-\Ml v) ) 



K) ~ 1 



In 



A 2 -t\l-T 



(37) 
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L(t) = 



In 



A 



2 



<2o 



(38) 



where the parameter A can be considered as a measure of the onset of the asymptotic 
behaviour. The value of Qq is strongly correlated to the one of A in the actual fits, we choose 
its value to be in the few hundred MeV region (as discussed in section [TV]), Q ~ Aq C d- We 
take into account three poles in the isovector (v) channels and denote these by p', p" and 
pin |g p con t r ib u tion j s fixed by extended unitarity as discussed in section |11 C| , 

**(* - 1 M> lmF ^ <L-\Ml) + blL~\M^{l-tKr^ 

t% [t) ~ ~ Lm$ L(t>)(t> -t)~ 2(1 - t/df) ' {M) 



71 



with the mass parameters M a and determined from a fit to Ff(t) and the af , 6f , cf , <if 
given after eq.(p3|). We find M 2 = 0.5 GeV 2 and M 2 = 0.4 GeV 2 . In the isoscalar channels, 
we have three meson poles, the first two with fixed masses at = 0.782 GeV and M$ = 1.02 
GeV plus one heavier state denoted by S'. The contribution of L(t) for momenta smaller 
than A is balanced by the appropriate normalization factors at the meson poles, L~ 1 (Mf), 

The analytic structure of the functions F?(t) in the complex t-plane is shown in Figs.3a,b. 
In the isovector case (Fig. 3b), I = (t>), the two-pion continuum leads to a cut from t = 4M 2 
to t up = 50M 2 . In addition, there are poles at t — M p r, M^', M p m which to some extent 
compensate the neglect of cuts related to higher mass continua like 4tt, NN and so on 
(with positive G-parity). At larger values of t along the real axis, there is a singularity at 
^sing = A 2 — Qq from the logarithm plus a cut from t sing to A 2 since the exponent 7 appearing 
in L(t), eq.(^), is rational. Finally, there is the right-hand cut from A 2 to 00. The only 
difference in the isoscalar case (Fig.3a), I = (s), is that one has three poles at M 2 , M| and 
Mgr. The latter two simulate cuts related to multi-pion and other intermediate states with 
negative G-parity. The singularity at t sing has no physical significance, it is related to the 
way in which the large — t behaviour is enforced. We will come back to this point in the next 
subsection. 

The normalization conditions of the form factors at t = take the form 

Pt ] S m + E ai = v{ L-\0) (40) 

with the v\ given in eq.© and |l 



1 /-50M2 I / \ 

- / 2 dt>t' h li*Fat')L-\t>) = -UL-\M*)(d^ (41) 

7T J AMI 2 \ / 



71 



* 4 This notation, however, should not imply a priori that these poles have to be identified with 
higher mass excitations of the p. This topic will be taken up again when we discuss the actual 
results of our fits. 
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The conditions on the neutron charge radius, eq. (|26|) , and on the proton charge radius, 
eq. fl29|) , translate into 



(o)i(o) = E 



a\L-\M})( 1 



+ 



7 



Ml \M] A 2 L-V7(0) 



7 



+ 2 



(42) 



If one only imposes the value of the neutron charge radius, one has to take the appropriate 



linear combination of the two conditions given in eq.(42). We stress again that this latter 
case should be considered as more realistic. The implementation of the pQCD constraints, 
eq . (|35D , via the logarithmic function L(t) is similar to the approach taken in ||T0| , pJ] , p6| . 
The following super convergence relations have to be fulfilled so that the proper pQCD 
asymptotics emerges 



m 



pf ) ^) + E«^" 1 ( M /) = o 



p ( 2 ) 5 IXv) + Y j m 2 i4l~\m 2 i ) = q . 



(43) 



These relations assure that the pole-terms and the 27r-continuum do only contribute to 
subleading orders as \t\ becomes large. In Appendix |A|, we discuss how all these constraints 
are technically implemented. 

Finally, let us count the number of free parameters for the fit functions eqs. (|36| , |37|) . For 
n s = 3 isoscalar and n v = 3 isovector mesons, one gets 2(n s + n v ) = 12 pole and (n s + n v ) = 6 
mass parameters plus the pQCD parameter A. We remark here that we do not treat Qq 
as a free parameter since only the combination L _1 (0) = [ln(A 2 /Qo)] 7 enters the pertinent 
formulae. We only make sure that Qq comes out in the few hundred MeV region because 
of its proportionality to Aqcd- We also stress that the p contribution does not induce any 
free parameter and that the mass of the lowest two isoscalar poles are fixed to M u and 
M^. With the 12 relations eqs.( ^0| , ^2| , ^3D the effective number of free parameters is thus 
18 + 1 — 2 — 12 = 5 (or 6 if we only enforce the condition on the neutron charge radius). 



B. Additive parametrization 



The form factor representation eqs. (|36| , |3~7D is very simple but one can not directly separate 
the hadronic from the pQCD contribution in the spectral functions. Therefore, we present 
here another form which is similar in spirit to the one proposed in ref. [pjj. Under some 
approximations, one can in fact construct a spectral representation for the pQCD part which 
can be added to the contributions from the meson poles and the two-pion continuum. This in 



12 



general introduces some additional fit parameters cf which are related to the normalization 
of the form factors at large t, 

J 

12M7 



~> 77^77^777^17 • (44) 



t*+i [hi(-t/Q 

In addition, one has to take special care about the singularity at t S i ng as detailed in ap- 
pendix [B]. It leads to an extra pole term not considered in |L7j |Tj|. The explicit additive 
parametrization takes the form 

~QCD , Tm pQCD/,,^ 



Im^ QCD (i) 



A 2 - Qo - t vr 7a2 i' - i 

/ r 7r 
sin 7 arctanf— — — ) + 7r0(l — 2) 
V L 111(2) 



(45) 



with z = I (A 2 — t) / Qq I ■ Clearly, the representation for the Im D is more complicated than 
the one used in |H] [0 but it leads to the desired behaviour, eq.(pH|). The normalization 
and radius constraints as well as the super convergence relations are modified, 



1 /,QCD 



V m 2 A 2 - q; 

(0) r , J 1 /,QCD . I t (0) n 



P? 5i M + E m? 4 + (a 2 - Ql) 4- QCD + 4 $ = , 

I,QCD r ,-pl 

+ -;d- 2| + E^ + ^»pS- 2, = ^(o) , (46) 



(A 2 - Q 2 ,) 2 1 lsl Y M ; 4 ' A 



with 



£<* > = i jT ^' f fe Imif CD (*') . (47) 

Here, the pf° are defined as in eq.(f4"T|) with L(t) = 1, and the a[' QCD are defined in ap- 
pendix |B[ The pertinent numerical values are dFi S \o)/dt = 1.37 GeV -2 and dF± v \o)/dt = 
1.30 GeV -2 . These last two conditions in eq. ffilf ) have to be combined appropriately if only 
the neutron charge radius as measured in low-energy neutron-atom scattering is enforced 
(i.e. only one constraint results). In the additive parametrization, eq . fl45|) , we have in prin- 
ciple threep 5 ! more effective free parameters than in the multiplicative one. In the latter 



* 5 Since we take the value of A 2 as given from the multiplicative parametrization, only three of 
the c{ are independent. 
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case, the normalization of the form factors at large momentum transfer are essentially fixed 
by our choice of the functions L(t). We stress here that the values one obtains for the pa- 
rameters c\ are only indicative of the strength of the form factors in the asymptotic region 
since their numerical values are very sensitive to the number of meson poles at low energies 
one accounts for (see also ][L7| |0|). In fact, what can happen in the additive parametriza- 
tion is the following. If one does not give reasonable bounds on the new parameters cf, 
a simple ^-minimization can lead to unphysical results with very large \c\\. These will 
then influence the low-i behaviour of the various form factors, in particular the isoscalar 
radii. Therefore, we constrain the additive parametrization to essentially give the same low 
momentum description of the four em form factors as does the multiplicative one. This also 
means that the c\ are fixed, i.e. they are no longer free parameters (see also app.[FJ). This 
method ensures that we can make sensible statements about the pQCD contributions to the 
various spectral functions. Only if one would have data at larger Q 2 it would make sense 
leaving the c\ as free parameters. 



C. Fits with an effective p pole 



The inclusion of the p contribution as detailed in section [II C| leads to a large value for 
the tensor-vector coupling ratio (as defined in eq. ([l5|) ) fl22f , 

k p = 6.6 ± 1 . (48) 

This value is in agreement with other determinations, see e.g. Grein [II] (n p = 6.1 ± 0.6). 
More recently, Brown and Machleidt [4H] have discussed the evidence for a strong pNN 
coupling from the measurements of the €\ parameter in NN scattering. In ref. [17|] an 



effective p pole with a mass of 0.63 GeV was used and led to typical values of k p = 5.9. We 
will also perform such a simple pole fit, i.e. substituting the full two-pion continuum by 
a p pole with a variable mass and taking in addition two more pole terms in the isovector 
channels. Our motivation to perform such types of simplified fits is to check whether the 
large value of the tensor-to-vector coupling of the p with the correct implementation of the 
pQCD constraints can be considered a generic result. 



IV. RESULTS AND DISCUSSION 

Before discussing the specific results of our fits, we wish to make some general comments. 
We had argued before that the masses of the three isovector excitations and of the highest 
isoscalar one need not to coincide with masses of physical particles. However, we have found 
that it is possible to find a minimum in the x 2 hyper-surface fixing M p > = 1.45 GeV and 
Mp/i = 1.69 GeV together with M$> = 1.60 GeV. These are the values of the most recent 
particle data group compilation [[£| for the lowest isovector-vector and isoscalar-vector 
meson excitations. Leaving the values of these masses completely free does not alter the 
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X 2 significantly. In contrast, the mass of the third pole in the isovector channel is tightly 
bound due to the various constraints the fits have to obey. We observe that the mass of 
this third isovector pole tends to come out close to one of the other isovector poles, thus an 
effective double-pole structure around Mp> or M p » emerges. This is in marked contrast to 
the findings of ref . |2l| . However, it is mandatory to retain three poles besides the two-pion 



cut contribution in the isovector channel. Also, independent of the details of the fits, we 
find that while the form factors Fi S \t) and F% < (t) exhibit a stable dipole structure (i.e. the 
lowest two poles have residua which are equal in magntiude but with different signs), this is 
not the case for F 2 {s \t) and F} v \t). These findings agree with the ones of [pi)] . Concerning 
the accuracy of our fits, all normalization, radius and superconvergence relations are fulfilled 
within machine accuracy, typically much better than one part in 10 12 . After these general 
remarks, we turn to a more detailed description of our results. 

A. The best fit: Form factors, radii and coupling constants 

The optimal fit to the available set of form factor data is obtained with the isovector 
masses of 1.45, 1.65 and 1.69 GeV, respectively and the isoscalar ones being 0.782, 1.019 and 
1.60 GeV. The corresponding residua are a± = —3.465, a 2 = —6.552, a± = 40.26, a 2 = 
7.881, a(" = -37.30, af = -2.821, af = 0.747, a% = -0.122, af = -0.738, af = 0.162, 
af = -0.0382 and af = -0.0406. We find A 2 = 9.73 GeV 2 and Q 2 = 0.35 GeV 2 . We 
notice that to have good fits, we can vary A 2 between 5 and 15 GeV 2 without drastically 
changing any of our conclusions. We always work with 7 = 2.148 since the fits are completely 
insensitive to the possible variation in this quantity. In this fit, only the constraint on the 
neutron charge radius is imposed and the p — uj mixing in the two-pion spectral function is 
included. The resulting x 2 /datum is 1.09. 

In Fig. 4, we show the electric and magnetic proton and neutron form factors normalised 
to the dipole fit.|^| Similar to earlier findings, we note that there are substantial deviations 
from the dipole fit in all channels, particularly at large momentum transfer. We also note 
that a better data basis is clearly needed. 

Of particular importance is the determination of the nucleon radii. In table 1, we give 
radii corrsponding to the Pauli and Dirac form factors in comparison to the results of 
ref. For the isovector form factors, the radii are indeed dominated by the two-pion 



plus p contribution, we have r{ ~ r 2 ~ 0.75 fm. The corresponding neutron and proton 



radii are given in table 2. The uncertainty for these radii is 1% (for comparison, in |2T 
the uncertainties on the radii were of the order of 3%). This number is calculated in the 
following way. In the parameter-space we look for solutions with a comparable x 2 /datum 
than the best fit has. Equivalently, one can sum in quadrature the la deviations of these 



#6 In the case of we divide by as given by Platchkov et al. [Q| (denoted by Gp) since in 
the dipole approximation, G% = 0. 
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parameters contributing to the various radii. Our results are comparable to the ones of ref. 
pl| with the exception of and r^. This can be traced back to the fact that in |2T| the 
superconvergence relation 

\a^f + £ 4il4 = 0, (49) 

/ V=p',p",p"' 

was not taken into account. In contrast, our values for and are based on a completely 
consistent calculation. We point out that there exist some on-going activity e.g. at ELSA 
(Bonn) to determine the neutron magnetic form factor more precisely at low and moderate 
momentum transfer. We also note that the value for r E is on the low side of the result 
of ref. |15||. If one insists on their central value, r E = 0.862 fm, by imposing the proper 



slope condition, one is not able to simultaneously fit the neutron and the proton form factor 
data. In fact, in |y| only proton data were considered. We conclude that the uncertainty 
attributed to r E in |15| is presumably underestimated. As we already anticipated in section 



! IE], the constraint from the proton charge radius in its present form is too restrictive to be 
applied in the dispersive analysis. Better low-energy data are clearly called for to settle this 
issue. 

The mesonic coupling constants can be directly inferred from the pertinent residua. 
Applying the same analysis to get a handle at the uncertainties as described above, we find 

g u,NN = 2QM ± Q 25 ^ g uNN = _ 3 A1 ± Q 2A ^ 

g *NN = _g 16 ± Q 23 ^ g *NN = 2 Q1 ± Q 33 ? (5Q) 



which are compared with the findings of ref. [21] in table 3. Our ^-couplings are somewhat 



larger but consistent within error bars. We note that the uNN coupling is larger than in 



typical one-boson exchange potentials or from the analysis of Grein and Kroll [j|5] using 
forward dispersion relations for AA-scattering, (gf NN ) 2 /4iT = 8.1 ± 1.5. Such a small 
coupling constant value cannot be accomodated in our fit, it is inconsistent with empirical 

(s) 

information on the slope of F x if the u and the $ lead to the dipole structure as described 



above. This point is also discussed in some detail in [27J. Furthermore, we remark that 
a direct comparison with coupling constants obtained in boson-exchange models, which in 
general include strong meson-nucleon form factors, has to be taken cum grano salis-P 7 !- In 
flavor SU(3), one can derive the following formulae for the u — $-mixing angle, O, 



V3 gf N gf NN 



NN - tall ( ) = -^JNN ■ ( 51 ) 



cos(Q)gf NN v ' g$ 



* 7 The photon couplings through vector mesons to the nucleon need not be the same than the 
purely strong interaction VNN couplings 
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Using g{ = 2.6 [p2], we have = 35°, very close to the ideal mixing angle of 37°. This 
means that the $ is almost entirely an ss state and is thus supposed to decouple from the 
nucleon (to leading order in flavor perturbation theory). This is the much discussed violation 
of the OZI rule. This apparent paradox finds its resolution in the fact that the simple pole 
approach for the QNN coupling effectively includes contributions from the KK continuum. 
Stated differently, the $ can couple to the kaon cloud surrounding the nucleon (as modeled 
e.g. in [p6fl). This topic could be investigated further along the lines discussed in section 



11 □ , i.e. by analyzing the NN — > KK partial waves. 

We close the discussion about the coupling constants with some remarks on n p . As argued 



in section [111 Q one can also perform fits with a p-pole. In this case, one cannot take the 



physical mass for the p since otherwise the isovector radii are severely underestimated. We 
have performed such fits and find 

re, = 6.1 ±0.2, (52) 



which is consistent with previous determinations as discussed in section [III Q We consider 



this an important consistency check on our fits. The realistic fits, however, have to include 



the correlated two-pion exchange as described in section |I C . 

Next, we discuss the large momentum behaviour of the form factors. Only for G P M , 
there are data for Q 2 > 10 GeV 2 . In Fig. 5, we show the quantities Q 4 G P M (Q 2 )/ fi p and 

L~\Q 2 )Q*G P M(Q' 2 )/(jh L ~ X ( )) U P t0 Q 2 = 50 GeV2 - Within the uncertainties, the curve 
representing the second function tends to a constant (as it is expected from pQCD), but 
it is obvious that more precise data at high momentum transfer are mandatory to really 
pin down this behaviour. The available data do not exclude that asymptotia sets in much 
beyond Q 2 = 30 GeV 2 . Often considered is also the ratio Q 2 F 2 (Q 2 ) / F^Q 2 ) which should 
become constant as Q 2 becomes large (the extra Q 2 in front of F 2 compensates the spin-flip 
suppression ~ 1/Q 2 ). In this ratio, many uncertainties related to the exact form of the 
nucleon wave function drop out. As seen in Fig. 6, the presently available data are at too low 
momentum transfer to test this prediction, although there is some hint of Q 2 ^(Q 2 ) / Fi(Q 2 ) 
becoming constant as Q 2 increases. Again, more accurate data at higher Q 2 are called for. 



B. Additive parametrization: Hadron versus quark contributions 

As explained before, the multiplicative parametrization is not well suited for separating 
the hadronic (pole) from the quark (pQCD) contributions. That is the reason underlying the 
additive parametrization. In this case, the four normalization constants q (of which only 
three are independent) in principle increase the number of tuneable parameters. However, 
as explained before, one has to impose certain constraints on the actual values of the q since 
otherwise completely unphysical solutions of the fitting procedure can emerge. Our strategy 
is therefore to constrain pole parameters and the c\ such that we essentially recover the 
low momentum description of the multiplicative parametrization, in particular the nucleon 
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radii. Only with such constraints one can make sensible statements about the separation of 
hadronic and pQCD contributions. 

In Fig. 7, we show a typical result for the isoscalar and isovector form factors. We have 
= -53.40, c[ s) = -8.20, 4 v) = 95.56 and cj> s) = -1.55. The x 2 /datum is 1.86. This 
increased value is a mostly reflection of the the approximations performed (i.e. the neglect 
of subleading 1/t corrections) to derive the additive parametrization, compare app.pj.p 8 ! We 
repeat that leaving all parameters free, one could naturally find a solution with a lower x 2 
than for the multiplicative parametrizations. Such solutions, however, have to be discarded 
as discusssed before. For all form factors, the hadronic and the quark contribution are 
of comparable size (in magnitude) around A 2 ~ 10 GeV 2 . We notice that for the Dirac 
form factors Fi S ' v \Q 2 ) } the hadronic contribution quickly drops off beyond Q 2 larger than 
10 GeV 2 . For the Pauli form factors F% f (Q 2 ), this fall-off is slower which essentially is 
the reason that one does not observe pQCD scaling for the available data. In all cases, 
the quark contribution is very small at low t (by construction). A similar behaviour was 
noticed in |17]], [18|| although in these papers the asymptotic behaviour was incorrectly 



implemented in the spectral functions. Of course, with the presently available data base, 
these results should only be considered indicative. With better data in the few and many 
GeV region, one will eventually be able to more cleanly separate the hadronic from the 
quark contribution. In particular, for the range of momentum transfer available at CEBAF, 
one will essentially probe the transition region from the hadronic to the quark description. 
The planned experiments at CEBAF [[|IJ will certainly shed light on this interesting regime. 
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* 8 In that appendix, it is also shown how one can systematically improve this procedure. For our 
purpose, the lowest order approximation used here is, however, sufficient. 
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APPENDIX A: IMPLEMENTATION OF THE CONSTRAINTS 



In this appendix we show how to evaluate the fit constraints. In general we have a system 
of 12 equations. It can be written in the following way, 



MlA\ 



CI 



(Al) 



where the matrix M/ has the form 



M 



(Mf 2 Mf 2 Mf 2 

1 1 1 
V M h M h M h 



Mr 



a ■ 



i 

Mf 

M 2 



7 L 1 /7(o) 



(A2) 



Setting a = 1, 0, respectively, permits to switch off certain vector meson pole terms, k 
specifies the parametrization. n = 1 if the multiplicative parametrization is used and k = 
for the additive one. With A\ = a I i L~ 1 (Mj) we can write 



Al 




CI 



p) 



\ 



(o)x 

& 



J 



(A3) 



— * j 

The vector C- contains all fit constraints. The first component enforces the correct normal- 
izations, the second and the last one in the case (i = 2) induce the superconvergence relation 
(fj3|). For (i = 1) the third component implies the slope informations stemming from the 
experimental values of the electric radii of proton and neutron, 



C 



where 3? is the contribution of the two pion continuum to the radius constraints 



(A4) 



ft 



7L 1/7 (0)\ XL- 1 {Ml 



2A 2 



(A5) 



We simply invert the matrix M\ to find the solution vector Af, 



{Ml 



1 

m 



( M h - M h Mf 3 2 M h - Mj 2 M h Mf 2 2 
M h - M h Mj 2 M h - Mj 3 2 M h My 2 
V M h - M h Mj 2 M h - Mf 2 M h MJ 2 - Mj* 



M T 2 
MT 2 



(A6) 



with 



m 



(Mf 2 - Mj 2 )M h + (Mr 2 - Mj 2 )M l2 + (Mr 2 - Mj 2 )M J 



(A7) 
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APPENDIX B: DERIVATION OF THE ADDITIVE PARAMETRIZATION 



Here, we derive the additive parametrization, eq.(^). The Cauchy integral representa- 
tion for the fit functions Fl(t) (i = 1, 2; I = (v), (s)), eqs. (|3^j37| ), takes the form 



dt 



2m Jci V - t 



(Bl) 



with C/ a closed integration contour as shown in Figs.3a,b. Cj is chosen such that Re F*(t) 
can be calculated as a Hilbert-transform. For that, we need the imaginary part of F?(t), 



Im Flit ± it) = Re F*{t) Im L(t ± z'e) + Im F*(t ± z'e) Re L{t) . 



(B2) 



The real and imaginary parts of F*(t) for t < t sing = A 2 — Q$ and t > t sing can be deduced 
easily. L(t) is purely real for t < t S mg- Consider now the imaginary part of L(t). For 
^sing < t < A 2 , we find (on the first Riemann sheet) 



Im Lit ± ie) = =F 



sin(77r) 



with 



A 2 



Ql 



(B3) 



(B4) 



and for t > A 2 , it reads (on the first sheet) 



sin 7 



Im L(t ± z'e) = =F- 



arctan( 



7T 



ln(z) 



+ 7T0(1 - Z) 



and finally 



(ln(z)) 2 + 7r 2 



Im L(A 2 ± ie) = , 



7/2 



(B5) 



(B6) 



assuming that Mf < A 2 for all /. Due to the singularity at t sing , special care has to be taken 
of the region t E [A 2 - 2Ql, A 2 ]. We have 



r / n 1 f 00 Re F/(t') Im L(t) 1 
Re if (t) = - / dt' — + — Im 

' w 7r Aa f - t 2tt 



f i v 



losing 
U 7 



+ ... , 



(B7) 



where the ellipsis stands for the contributions from the pole terms and the two-pion contin- 
uum (which are of no relevance for the large-t behaviour discussed here). For the contour 
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C" ng we choose a closed circle around the singularity with radius Qq. In the region of inte- 
gration and close to the singularity, we have t 3> Mf so that we can approximate Ff(t) by 
its asymptotic behaviour ™. 



Re Flit) = ± 



dt 1 



Im L(t) 



+ — Im 



7T 7 A 2 (t'Y+ 1 (t - t) 2TT 



dt' 



L(t') 



+ . . . . 



(B8) 



Let us concentrate on the second term in eq. (|B8|) . Parametrizing the integration path via 
t(<f>) = t sing + Qq exp(i(f)), we find 



I(t) 



1 r 2n 



iQge**ln(l -e 



2tt Jo " r [A 2 - Ql{\ - e^)f +1 [A 2 - t - Q 2 Q (1 - e<*)] 
We expand the first term in the numerator in eq.(B9|), 

i i / Ql 



(B9) 



[A 2 - Q 2 (l - e^)] i+1 [A 2 - Ql] l+1 V y " ' " 7 A 2 - Q 2 
and drop the term proprotional to exp(i0) in the numerator of eq. (|B9|) . Altogether, we find 



1 



1 



e i0 + . 



(BIO) 



I(t) 



~QCD 



a 



A 2 - Q 2 - t ' 



(Bll) 



with 



~QCD 



Ql 



2tt(A 2 - q 2 ) i+1 



C {1) - (i + 1) 



A 2 -g§ 



c( 2 ) + . . . 



2- 



dc/>e m ^ ln(l -e # ) 



<M-7 



(B12) 



where the coefficients C^ n ' can be found by numerical integration. They are purely real. It 
is sufficient to keep the first two terms in the series with = —5.861 and = 7.129. 
Setting now 



/,qcd / ~QCD 



(B13) 

we have arrived at the desired result, eq. (|45l) . 

Finally, we wish to establish a more rigorous derivation which is useful if one wants to 
retain more terms of the asymptotic expansion of the F?(t). For that, we expand the F*(t) 
in powers of 1/t, 



F!(t)= £ 4r ki 

fei<(i+l) 



(B14) 



#9 In 

our fits, the most massive isovector (isoscalar) pole is at 2.8 (2.6) GeV 2 , which is a bit close 
to the value of A 2 ~ 10 GeV 2 . This induces some uncertainty due to the terms neglected. 
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with the coefficients c{. given by 



ki 



j 7T J4M$ 



and 



c 



fci<(i+l) 







because of the super convergence relations. This leads to 

' 1 r°° . , Im Lit) 1 

1 — Im 



fci<(i+i) 



(t') fcl (# - t) 2tt 



dt 



, L{t>) 



c s ing (^fcj (J, _ t ) 



which can be written in a form similar to eq. 



ReFi(t) 



/,QCD 



+ - / dt 



A 2 -Q 2 -t ' n J A 2 

but in this case we find 

oo 



t'-t 



+ ... 



2tt(A 2 - Ql) k > 



Ql 



(B15) 
(B16) 



(B17) 
(B18) 



fci<(i+l) 

Imi^ QCD (t) = ReF^t) ImL(t) . (B19) 

Considering only the first term of the sum over ki in eq. (|B19|) is equivalent to the approxi- 
mation eq . ( p8|) . 



APPENDIX C: PARAMETRIZATION OF THE BEST FIT 



Here we give the parametrization of the best fit discussed in section |1V A| explicitly for 
easy usage, with Q 2 in GeV 2 , 
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Tables 





7i [fm] 


[fm] 


ri [fm] 


T2 [fm] 


Best Fit 
Ref. |21| 


0.782 
0.77 


0.845 
0.837 


0.765 
0.76 


0.893 
0.863 



Table 1: Radii of the Dirac and Pauli form factors. 





r p E [fm] 


r p M [fm] 


r\. L [fm] 


r\ [fm] 


rf [fm] 


[fm] 


Best Fit 
Ref. 21 


0.847 
0.836 


0.836 
0.843 


0.889 
0.840 


0.774 
0.761 


0.894 
0.883 


0.893 
0.876 



Table 2: Proton and neutron radii. 





{gf NN f/^ 




(at™) 2 /** 




Best Fit 
Ref. 1 21 1 


34.6 ±0.8 
30 ±3 


-0.16 ±0.01 
-0.17 


6.7 ±0.3 
4.4 ±1 


-0.22 ±0.01 
-0.3 



Table 3: Coupling constants of the isoscalar vector mesons. 
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Figure captions 



Fig.l Two-pion cut contribution to the isovector nucleon form factors. 



Fig. 2 Spectral distribution of the isovector form factors weighted with 1/t 2 . Shown are 
2ImF?(t)/t 2 (long-dashed line), 2rImF 2 " (t) /t 2 (short-dashed line) and 2lmG v E (t) /t 2 
(solid line). Upper panel: no p-u mixing. Lower panel: with p-u mixing. 



Fig. 3 Analytic structure of the fit functions eqs.( |36| , |3"TD in the complex t plane, (a) Isoscalar 
case. (1) Poles at M u , M$ and M S /, (2) the singularity at A 2 — Ql and (3) the right- 
hand cut starting at A 2 . The data are given at negative t. (b) Isovector case. The 
right-hand cut (1) is the two-pion continuum, (2) are the three isovector poles, (3) 
and (4) are equivalent to (2) and (3) of the isoscalar case. 

Fig. 4 Optimal fit within the multiplicative parametrization including p-u mixing and the 
constraint from the neutron charge radius. The data for G% are from |5(J |^7| [[KJ |49j 
[0 g @ H, for G P E from for G P M from 





n m m mi |9 



and for G n M from @ g @ @ ||7|. For the 
data of ref. [43], we have taken the values based on the Paris potential but enlarged 
the error bars to account for the model-dependence. 

Fig. 5 Asymptotic behaviour of the proton magnetic form factor. The best fit in com- 
parison to the available data is shown. Upper panel: Q 4 Gf / (Q 2 )//i p , lower panel: 
L~ 1 (Q 2 )Q 4 G^ [ (Q 2 )/(p p L~ 1 (0)) 

Fig.6 The ratio Q 2 F|(Q 2 )/Ff(Q 2 ) for the best fit compared to the data. 

Fig. 7 Hadronic (short-dashed lines) versus quark (long-dashed lines) contributions for the 
isoscalar and isovector Dirac and Pauli form factors (solid lines) in the additive 
parametrization. The absolute values of the various contributions are shown. 
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Figure 3 
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Figure 4 
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Figure 6 
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